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Abstract 

We study non-standard top quark couplings in the effective field theory approach. 
All nine dimension-six operators that generate anomalous couplings between the elec- 
troweak gauge bosons and the third-generation quarks are included. We calculate their 
contributions at tree level and one loop to all major precision electroweak observables. 
The calculations are compared with data to obtain constraints on eight of these oper- 
ators. 



1 Introduction 



Top quark interactions could provide relevant information on physics beyond the Standard 
Model (SM). Anomalous top quark interactions at colliders have been studied in the literature 
(see, e.g. Pfl|2lElSll5lEllZlElElliniinilI2llI3] and references therein). In particular, a 
model-independent approach based on a low-energy effective field theory is used to describe 
possible new physics effects. In this approach, high-scale physics is integrated out to obtain 
effective interactions that involve only the SM particles. These interactions are suppressed 
by inverse powers of A, the scale at which the new physics resides. 

Such a field theory must satisfy the SU(3)c x SU{2) L x U(l)y symmetry of the SM. 
With this requirement, the only possible dimension-five operator violates lepton number 
conservation and is irrelevant to top quark physics [14J . Thus the leading effects are generated 
by operators of dimension-six: 

C eS = C SM + ^ ^(C* + h.c), (1) 

i 

where Oj are the dimension-six operators and Cj are dimensionless coefficients. A complete 
list of dimension-six operators was given in [15j [16j [17]. Subsequently it was found that 
several of this operators are not independent [HUB]. A list of 59 independent dimension-six 
operators is given in [19] . 

New physics that affects top-quark interactions can be parametrized using these oper- 
ators. In particular, consider anomalous top quark couplings, such as the Wib, Ztt and 
'yit vertices. There are nine dimension-six operators that can modify the couplings of the 
third-generation quarks to the W, Z and 7 bosons: 
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where q is the third-generation left-handed quark doublet, t and b are the right-handed top 
and bottom, <fi is the Higgs boson doublet, e = ( ^ ) ' ^ = 6< ^*' an< ^ ^ M = ^" ~~ 

i^T T W^ — ig'YB^ is the covariant derivative where r l denote the Pauli matrices. WL, = 
— d v W* + geuxW^W^ and B^ u = d^Bl — duB 1 are the field strength tensors of the 
W and B field. The contribution of these operators to the vertices can be found in [18] . 

Naively, from dimensional analysis we may expect that the effects of these operators 
are suppressed by E 2 /A 2 where E is the energy scale of the process. This is not the case 
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for the operators listed above in Eqs. (|2|)-( TT0|) . The anomalous couplings generated by these 
operators violate the SU (2) L symmetry, so they are related to the Higgs vacuum expectation 
value v. Instead of E 2 /A 2 , these anomalous vertices scale as v 2 /A 2 , which is independent 
of the energy scale of the process. This can be seen in [18], where the relation between the 
anomalous couplings and the dimension-six operators are given. 

The consequence of this scaling is that the effects of new physics will not increase with 
energy. On the other hand, the effects will not disappear in the low energy limit. Therefore 
an important question is whether it is possible to extract better bounds from electroweak 
precision measurements for these operators, than from the measurements performed at high- 
energy colliders. 

The electroweak precision measurements have a much cleaner background than hadron 
colliders, and therefore are performed with a higher level of precision. However, most of the 
operators listed above do not directly contribute to these measurements at tree-level. Their 
corrections to the W, Z and 7 self-energies occur at loop-level, so they are suppressed by 
g 2 j (47r) 2 . Still, the large mass of the top quark can lead to an enhancement of the loop-level 
contribution, and as a result the constraints on top-quark anomalous couplings obtained from 
precision measurements may be comparable with those obtained from collider experiments. 
The top quark plays an important role as a virtual particle in precision electroweak physics. 
Indeed, the correct range for the top-quark mass was anticipated by precision electroweak 
studies. Now that the top-quark mass is accurately known from direct measurements, we 
can ask what the precision electroweak measurements have to say about the presence of 
dimension-six operators in loop diagrams involving the top quark. 

The easiest way to put constraints on these operators is to consider the "oblique pa- 
rameters" [2"CT| |2"T| |2"2"] . as most of the operators contribute only through corrections to the 
self-energies of the electroweak gauge bosons. However, as we will see in Section 2, this ap- 
proach is not appropriate for all nine operators. Therefore we explicitly calculate the effects 
of these operators on all electroweak measurements. We compare the results with data and 
perform a global fit to obtain one-sigma bounds on the coefficients of these operators. In 
order to be specific, we assume that these operators (plus two additional operators) are the 
only new physics effects in the theory. 

The operators in Eqs. (l2"j)-f lT0|) may be expressed in any basis that is convenient. We 
choose a basis in which all fermion fields are mass eigenstates, with q = V^l + VtsSL + 
Vtdd>L) [23]. In W boson self-energy diagrams, one must sum over all charge -1/3 quarks. 
These quarks are much lighter than the top quark in the loop, so it is a good approximation 
to neglect their masses, in which case the CKM factors add up to unity. The only place 
that a CKM factor survives is in photon and Z boson self energies involving a 6-quark loop 
and the operators Obw or ObB- These diagrams yield a factor of V t b, which is very close to 
unity for three generations, and can be ignored. The operators in Eqs. fT2"])- (fT0]) also give 
rise to nonstandard effects such as flavor- changing neutral currents, right-handed charged 
currents, etc. [21]. We do not explore the constraints on the operators from bounds on these 
processes. 

The paper is organized as follows. In Section 2, we discuss the oblique parameters and use 
them to constrain one operator as an illustration. In Section 3, we show all major precision 
electroweak measurements that we will use to obtain bounds. In section 4, we calculate the 
corrections to all observables from these operators, and perform the global fit. We present 
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our conclusions in Section 5. Finally, we show the self-energy corrections from each operator 
in Appendix |Aj and some numerical results of the global fit in Appendix Q31 

2 Oblique parameters 

The operators listed in Eqs. fl2|)- f|T0|) affect the precision electroweak measurements in two 
different ways: 

• > and O^, modify the Z — >■ 66 measurements at LEP at tree-level. 

• All operators modify the self-energies of W, Z and 7 at loop-level, and therefore affect 
all measurements indirectly. 

The first effect is equivalent to a correction to the Zbb couplings: 

^ = -~(cS' + C«), (11) 

s& = -It* **' (12) 

where g b L , g b R are the left- and right-handed Zbb couplings. This correction is easily included 
in the calculation. 

Operators listed in Eqs. ff2|)- f|T0l) modify the gauge boson self-energies through the loop 
diagrams shown in Figure [flj They affect all quarks and leptons universally through gauge 
boson self-energies. This kind of effect is referred to as "oblique". Traditionally three 
parameters, S, T and U, are used to describe the oblique new physics [201ET]. The idea is to 
Taylor-expand the four self-energies 11^/^) ^-zz, n 77 and IL^, which only include the new 
physics contributions, to order q 2 . Requiring the photon to be massless, Il 77 and U^z must 
be zero at q 2 = 0, so there will be six non-zero coefficients. Three of them are absorbed in 
the definition of g, g' and v. This leaves three independent parameters. The definitions of 
the S, T and U parameters are given in |25j . 

We list the contribution of the nine dimension-six operators to the gauge boson self- 
energies in Appendix [A] The coefficients C{ may be taken to be real because an imaginary 
part of Ci violates CP and will not contribute to any self-energy. It is straightforward to 
calculate the three oblique parameters using these expressions, and to compare with the 
existing bounds on the S, T and U parameters. 

We presented this analysis for the operator O t w in Ref. [26]. We found that the S 
parameter is divergent, which is not surprising because in an effective theory there are two 
dimension-six operators that contribute to the S and T parameters, respectively, at tree 
level: 

WB = (0VVX,^, (13) 
Of = (0tl>V)[(iV)V]- (14) 

1 There is also a diagram contributing to the VK-boson self energy, with a top-quark loop, constructed 
from the four-point contact interaction given by O t w and O^w- Since this interaction is antisymmetric in 
li, v, this diagram does not contribute to the self energy. 
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Figure 1: Corrections to gauge boson self-energy. The black dots indicate the dimension-six 
vertex. 

The divergent terms in the S and T parameters are absorbed by renormalization of the 
coefficients of these two operators, so S and T do not give useful information on the size 
of On the other hand, the U parameter is finite, because there is no dimension-six 

operator that contributes at tree level. Using the experimental bound U = 0.06 ± 0.10 given 
in [25], we found 




Unfortunately, this kind of analysis is not appropriate for all nine operators. The S, 
T, and U parameters are defined by assuming a linear q 2 dependence of the self-energies. 
However, once loop-level contributions are included, the self-energies contain terms like In q 2 
and q 2 lnq 2 . In particular, in a diagram with a bottom quark loop, the self-energies can 
have very different q 2 dependence in the regions q 2 < Am 2 and q 2 > Am 2 . An example is 
shown in Figure EJ Since the precision electroweak measurements include data measured 
at both q 2 rs and q 2 > m|, it is not reasonable to use a bound obtained by assuming a 
linear q 2 dependence. In addition, this calculation does not make full use of the obtained 
q 2 dependence of the self-energies. By calculating the U parameter, one can only put a 
constraint on one special linear combination of the operators. The precision measurements, 
on the other hand, contain much more information. 

In order to fully study the effects of the nine operators, we will explicitly calculate the 
effect of self-energy corrections on all electroweak measurements, and perform a fit includ- 
ing all operators in Eqs. fTSl- fflOl and Eqs. f|T3l) - f|T4l) . The self-energy corrections given in 
Appendix [A] contain divergent terms. Since all divergent terms are either constant or pro- 
portional to q 2 , they can contribute at most through the three oblique parameters. Therefore 
all divergences can be properly absorbed once we include Owb and in our analysis. An 
effective field theory is renormalizable in this sense. 

2 We found that, rather than being divergent, the contribution of Otw to the T parameter vanishes [26] , 
A top-quark model that gives a nonvanishing contribution to the T parameter is discussed in Ref. 





0.7 ± 1.1 TeV~ 2 . 



(15) 
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Figure 2: The q 2 dependence of H Zr The contributions from the operators Of q and 0, b are 
shown for illustration. A linear part in q 2 is subtracted so that I1 7 ^(0) = n^ z (0) = 0. There 
is a branch point at q 2 = 4m 2 . 

3 Experiments 

The measurements we use to constrain the coefficients of the operators are listed in Table [TJ 
Detailed descriptions for individual experiments can be found in the corresponding references. 
For a given observable X, the prediction of the effective field theory can be written as 

Xt^XsM + ^^Xf 1116 , (16) 

i 

where X t h is the prediction in the presence of the operators, X$m is the Standard Model pre- 
diction, and Yli C%Xf im6 are the corrections from the new operators. Since only dimension-six 
operators are included, higher-order terms in Ci/A 2 are dropped. 

The SM predictions are computed to the required accuracy, and can be found in the 
literature shown in Tabled] The three most precisely measured electroweak observables, a, 
Gf, and mz, are taken to be the input parameters, from which the SM gauge couplings and 
the Higgs VEV are inferred. In addition, the following input parameters are used: 

m ffigg s = QOtll GeV > m t = 173.2 ± 1.3 GeV, a s (m z ) = 0.1183 ± 0.0015, (17) 

except for LEP2. The sensitivities of the SM predictions to the input parameters for the 
fermion pair production and W pair production cross sections at LEP2 are negligible com- 
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Notation 


Measurement 


Reference 


Z-pole 


T z 

Chad 

R f(f = e,H,T, b, c) 
ipn f f = e, u, t, 6, c, s) 

Af(/ = e,/x,r, &, c,s) 


Total Z width 
Hadronic cross section 
Ratios of decay rates 
Forward-backward asymmetries 
Hadronic charge asymmetry 
Polarized asymmetries 


mum 


Permion pair 
production at LEP2 


c/(/ = e) 
A f FB (f = »,T) 


Total cross sections for e + e~ — > ff 
Forward-backward asymmetries for e + e~ — > // 


ED 


¥ V A A A CAO O 

and decay rate 


111 \y 

r w" 


W mass from LEP and Tevatron 
W width from Tevatron 


25 


DIS 

and 

atomic parity violation 


Qw(Cs) 

Qw{Tl) 
Qw(e) 
9l,9r 

g v v e ,9A 


Weak charge in Cs 
Weak charge in Tl 
Weak charge of the electron 
i^-nucleon scattering from NuTeV 
v-e scattering from CHARM II 


m 


W pair production 




Total cross section for e + e~ — > W + W~ 


28, 



Table 1: Major precision electroweak measurements used in this analysis. The total cross 
section for e + e~ — > e + e~ is divergent. We use the cross section in the angular range cos# e 
[-0.9, 0.9] instead. 



pared to the experimental errors [29J. Therefore, we use the SM prediction given in the 
corresponding references. 

The corrections from the new operators include the tree-level contribution to the self- 
energies from Owb and the tree- level correction to the Zbb couplings from O^, 
and O^, and the loop- level contribution from all nine operators in Eqs. (l2"j)- ([Tu]) to the self- 
energies. Once the self-energies are given, the corrections Xf imG to all the experiments can 
be obtained from the modified tree-level formulae for each observable. This will be discussed 
in the next section. 

Given these results, we can calculate the total \ 2 as a function of Cf. 



X 



2 \ ^ {Xth — X exp ) 2 \ ^ (-^SM _ -^cxp + Yli CiXf lra& ) 2 



where X exp is the experimental value for observable X and ax is the total error which consists 
of both experimental and theoretical uncertainties. The \ 2 is a quadratic function of Cj. 
The fit for the coefficients of the new operators is given by minimizing \ 2 - The one-sigma 
bounds on the coefficients are given by \ 2 ~ Xmin = 1- 

Eq. (ITS]) needs to be modified to account for the correlations between different measure- 
ments. There are two sets of data for which the correlations between measurements cannot 
be neglected. These are the correlations between Z-pole observables [27], and the exper- 
imental error correlations for the hadronic total cross sections at LEP2 [2H]- To include 
correlations, Eq. (ITS]) should be modified to 

X 2 = £( X sm - X^ + J2 aX?' d ™ 6 )(a%\Xl M - X| xp + dim6 ) (19) 

p,q % i 

where X p,q denotes different observables. The error matrix a 2 is related to the error a p and 
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the correlation matrix p pq by 

°lq = VpPpqVq (20) 

The correlations for theoretical and experimental errors should be taken into account sepa- 
rately. 



4 Calculations 

In the presence of the new operators, the corrections to the self-energies of W, Z and 7 can 
be written as 

U XY = J2CiU XYl , (21) 

i 

where Hxy only includes the contributions from the new operators. (XY) = (Z Z) ,(WW) , 
(77), (jZ). 

For the operators in Eqs. fl2l- ffl0|) . the HxYi's are given in Appendix [S] We also include 
Owb and Oy in our calculation, so that the divergences can be absorbed. For these two 
operators, the contributions at tree-level are: 

IW = 0, (22) 

2v 2 v 2 
Rzz = C WB —s w c w q 2 + Cf ] —m 2 Z} (23) 

2v 2 2 

II 77 = — CwB-^;S W c w q , (24) 

v 2 

IT 7 z = -C WB — (c 2 w - s 2 w )q 2 , (25) 

where sw = sm9w and cw = cos 9w- 

In this section, we discuss the effect of self-energy corrections on each experiment. We 
will show how to obtain the CiXf ira% term in Eq. ffl8|) . We first illustrate the idea with an 
example. 

For processes involving light fermions as external particles, Peskin and Takeuchi have 
shown in Ref. [21] that the corrections to the gauge boson self-energies can be incorporated 
by a change in the couplings and gauge boson parameters. For example, for electromagnetic 
interactions, the coupling should be replaced by 

«,(g 2 ) = a (l + n; 7 (g 2 )), (26) 

where ao is the renormalized coupling, not including contributions from dimension-six op- 
erators. The self energy II 77 (g 2 ) contains only the contributions from the dimension-six 
operators, and II 77 (g 2 ) is defined as 



n 77 (g 2 ) = ^P-. (27) 



In the presence of a dimension-six operator, the renormalized coupling, «o, is different 
from the coupling measured in experiments. Therefore, the self-energy corrections affect 
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the theoretical predictions in two different ways, which we will call direct correction and 
indirect correction. The direct correction is simply described by Eq. (|26|) . Any observable in 
an electromagnetic process is affected by a change in the coupling. The indirect correction 
arises from the fact that we take the fine structure constant as one of the input parameters. 
The parameter ao is then shifted from the measured fine structure constant a, which is 
measured at q 2 = 0. Thus, from Eq. ( 126]) . 

a = a (l + n; 7 (0)). (28) 

Therefore any observable that depends on the fine structure constant as an input parameter 
is affected by Eq. fl28l) . We can now eliminate a by combining Eqs. fl26|) and ([28]) . to obtain 

a^q 2 ) = a[l + U^(q 2 ) - 11^(0)] , (29) 

which can be used to calculate the correction to any electromagnetic observable. 

We will show the direct correction and indirect correction to all observables in Section 
4.1 and 4.2, respectively, and combine them to calculate the total effects on all electroweak 
measurements, except for the cross section for W pair production. The W pair production 
cross section at LEP2 has relatively low statistics, and thus we will only consider the tree- 

(3) 

level contribution, i.e. the contribution from Owb and 0\ . 



4.1 Direct correction 

In the SM, the matrix elements of the charged- and neutral-current interactions mediated 
by electroweak gauge bosons can be written at tree level as 



GO 1 p 2 1 

-Mnc = e 2 ^ + „{h ~ s wQ)-2-—^2^ - s wQ')> (3°) 
q ^w^w " " l z 

e 2 . 1 



^ = ^7^'- (31) 



Peskin and Takeuchi have shown in Ref. [21J that the modification of the gauge boson 
self-energies can be included by writing 



.Mnc = e 2 M- + -T^H'a - 4^)^V( J 3 - *w*Q% (32) 
e 2 , Z v 



Mac = 2 /-, (33) 

2s 2 w , q 2 - m 2 w , 
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where the starred quantities are functions of q 2 : 





= (1 - Z w )q 2 + Z w (m 2 W0 + Ii ww (q 2 )) , 


(34) 


2 I 2\ 

m z AQ ) 


= (1 - Z z )q + Z z (m zo + Tl zz (q ) ) , 


(35) 


z w 


= 1 + ^ U ww{q 2 )\ q ^= m 2 w , 


(36) 


Z z 




(37) 


Zw*(q 2 ) 


= 1 + ^ww{q 2 )\ q ^ w - n; 7 (g 2 ) - ^n; z (g 2 ), 


(38) 


Z z *{q 2 ) 


= 1 + -^lW)| g2=m| - n; 7 (g 2 ) - ^^% z (q 2 ) 
dq z z " SwCw 


(39) 




= s 2 WQ - s w c w Tl^ z (q 2 ), 


(40) 


e 2 M 2 ) 


= e 2 + e 2 W^q 2 ), 


(41) 



where U' XY (q 2 ) is defined as 

W XY (q 2 ) = (U XY (q 2 ) - ll xy (0)) /q 2 . (42) 

The subscript denotes the renormalized parameter, not including contributions from dimension- 
six operators. When calculating the corrections due to dimension-six operators, we can use 
tree-level relations between the renormalized parameters, such as 



2 



2 _ e V 2 e Q V / A o\ 

m W0 — ~2 T> m ^0 - o ^2 7 ' 
b W0 ^ b W0 c W0 ^ 

Eqs. (132]) and (|33|) have exactly the same form as the tree-level SM amplitudes, except that 
all the couplings and gauge-boson parameters are replaced by starred parameters. This shows 
that the oblique corrections affect electroweak interaction observables only via the starred 
parameters. In other words, given an observable in terms of renormalized parameters at tree- 
level, we only need to replace the renormalized parameters with their starred counterparts 
evaluated at the appropriate momentum to incorporate the corrections from the self-energy 
diagrams. For example, at tree-level the left-right asymmetry A e at the Z-pole is given by 



AM) = 2(1 As \ * 



/' (44) 

l + (l-4^ ) 2 



This is modified to 



2(l-4,Um|)) 
l + {l-As 2 w M 2 z)Y 



after the self-energy corrections are included. Similarly, the Z to e + e partial width is now 
corrected to 

Note that these corrections come from the difference between quantities with subscript * 
and quantities with subscript 0, therefore these are direct corrections. 
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For low energy measurements, it is more convenient to write 

M NC = -±V2G F0 (l - -1-11^(0) ) (I 3 - s 2 w MQ) {I's ~ s 2 w MQ') , (47) 
V m z J 

Mcc = -2V2G F0 (l - -^-U ww (0) j (48) 

where 

G ™ = vb (49) 

is the Fermi constant, not including contributions from dimension-six operators. The direct 
corrections to any low energy observables are thus incorporated by replacing swo by sw*(0) 
and including an overall factor of (1 — H zz (0)/m z ) for neutral-current observables, and 
(1 — Hww(®)/ m w) f° r charged-current observables. 



4.2 Indirect correction 



The indirect corrections arise from the shifts in the renormalized parameters. The elec- 
troweak parameters (g, g', v) are not directly measured. Instead, we derive them from the 
most precisely measured observables (a, m z , Gp). When calculating the SM predictions for 
these observables, the tree-level relations between (g, g', v) and (a, m z , Gp) are used. When 
we include the new operators, the SM relations are altered. This corresponds to a correction 
to all three input parameters. 

To consider the indirect corrections, we use (cto, fn zo , GVo) to denote the renormal- 
ized electroweak parameters, not including contributions from dimension-six operators. The 
relation between a and «o can be read off from Eqs. ( l32i) and ( HIT) [see Eq. (128"]) ]: 



a 



e 2 M 
Ait 



a (i + n' (o)). 



The Z mass m z can be obtained by solving m Zit (m z ) = m|; this gives 



rrir 



m 



zo 



+ U zz {m 2 z ). 



The Fermi constant can be read off from Eq. ((48 

1 



Gp — G 



1 - 



-n 



m 



WW 



(0) 



w 



(50) 



(51) 



(52) 



Observables measured at energy scales above the Z pole are expressed in terms of (a, m z , 
sly) rather than (a, m z , Gp). To include indirect correction due to s^, we need 



"wo 




Kt (°) 



-n 



rn 



WW 



w 



(0) - -js-iwmi; 

1 1 L ry 



rrtr 



-U zz (m 2 z ] 



(53) 
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where t 

47ra 




V2G F m 2 z 



(54) 



is the value for calculated at tree level using the observed values for (a, m z , Gp). 

Combining Eqs. (|34l)- ()4Tl) with Eqs. (I50l) -(l53l) to eliminate ao, mzo, Gpo, and swo, we 
conclude that, for q 2 > 0, the total effect of direct and indirect corrections can be incorpo- 
rated by making the following replacement to the renormalized parameters in the tree-level 
expressions for any observable: 

{1 for interactions mediated by photon 

Zz*{q 2 ) for interactions mediated by Z bosor(55 
Zw*(q 2 ) for interactions mediated by W boson 

m% -> m%t = m% + 5m 2 z = m 2 z -Ii zz {m 2 z ) + Tl zz (q 2 ) - (q 2 - m 2 z )—Uzz(q 2 )\q2 =m % , (56 



'Ml S W* ~ S W $ S W ' *W 

s w 



c w b w \ " l w " L z / . 



The renormalized parameters (subscript 0) are now completely eliminated, and served only 
as intermediate quantities in the derivation. 

For any observable measured at the Z-pole or above, we can write it at tree level in terms 
of a, m| and s^: 

X ^ = X^(a,m 2 z ,s 2 w ). (58) 
Therefore the contribution from the self-energy corrections can be written as 

f) \^tree a \^tree a \^tree 

5X = CiX?** = ^^5a + ^Pt5m 2 z + ^^5s 2 w . (59) 
da dm z dsfy 

Note that Eq. ( 1581) is a tree-level relation, and we will not use it to compute the entire 
theoretical prediction. Instead, we use Eq. ( 1591) to find the corrections which arise from the 
new operators. Since these are already small corrections, the tree-level calculation is enough. 
We then add them to the SM predictions including radiative corrections, which are provided 
in the references shown in Table [TJ 

If the observables depend on the Zbb couplings, we will need to add to the r.h.s of Eq. (I59"|) 
the following terms: 



y2 f r (3) , r (l) , r \ dX lh° V 2 ( r (3) , Ml) r \ dXjjT / fin x 

- 2A2 [ C H + C H + - l 6 ^ + C ^ - 6 *>J ~d^- (60) 

This accounts for the tree-level correction to the Zbb couplings from Cf£j, C9^ and C^b, given 
in Eqs. (ITT]) and (fT2"]l . 

For low energy measurements, we can now write 

M NC = -Ay/2G F pM (J 3 - s 2 w MQ) (/a - SwMQ') , (61) 

Mcc = -2v / 2G F /+/_, (62) 
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where 

p*(0) = 1 - — U zz (0) + —Uww(0). (63) 

The results of DIS and atomic parity violation experiments are usually expressed in terms 
of the effective couplings in the neutral-current interactions. The corrections to these results 
can thus be obtained by replacing s 2 ^ by 



.(0) 



>w*V'i ■~ , \\ 



(64) 



and including an overall factor of p*(0) to the couplings. 
4.3 Observables 

Now we proceed to consider the correction to each observable. We will give the tree-level 
expressions for each observable, and then use Eq. ( 1591) to find the corrections that arise from 
the new operators. 

4.3.1 Z-pole process 

The process e + e~ — > f f was studied around the Z-pole at SLC and LEP1. At tree- level, 
the measured cross sections and asymmetries can be derived from two quantities: the partial 
width of Z — > ff, Tff, and the polarized asymmetry Af. The expressions are 

Af = -#4?, (66) 
9v +9a 

where the Z-fermion couplings g y and g* A are given by 



(67) 



The Z-pole observables include: 
• Total width 



/ 




A 


A 






+1 


+\ 


e,(i,r 


i 

2 


+ 2s 2 w 


l 

2 


u, c, t 




- ^s 2 


+ \ 


d, s, b 


i 

2 


+ ^s w 


1 
2 



Total hadronic cross-section 



r _ r ee rhad 

z 



< = (69) 
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• Ratios of decay rates 

f ^ for / = e, a, r 

for/ = 6 , c • (™) 

k r had J 

• Forward-backward asymmetries 

A FB = ^ A eAf, f = e,/i,T,b,c,s. (71) 

• Effective angle extracted from the hadronic charge asymmetry at LEP and from the 
combined lepton asymmetry from CDF and DO 

s 2 i = sV (72) 

• Polarized asymmetries 

Af, f = e,n,r,b,c,s. (73) 

With these tree-level expressions, we can apply Eq. (|59p to derive the correction from the 
new operators. For example, for the ratio Rb, we find 



M 164-36^ + 9 



(88^-84s^ + 45) 2 



fX z (m|) + (n; 7 (o) + -^n w (o) - -Ln^< 



V 



(3) , 40s w - 964, + 724, - 27 AQ v 2 203^ - 18s 2 w + 9 



24^ + " / T ' " " - 48^C, 



A 2 V * ' * J (88s 4 w - 8As 2 w + 45) 2 A 2 vb (88s 4 w - 8As 2 w + 45) 2 ' 

Using the expressions for the lily's given in Appendix^ this can be written in the form of 

In practice, instead of deriving the above equation for all observables, we actually did 
the following: for any given operator, we first give a small value to its coefficient Cj, then 
numerically compute all starred quantities using Eqs. fl55l 1561 1571) and fl63| |64|) . We then 
take the tree-level expressions for all observables and substitute in the starred quantities, 
to obtain the tree-level predictions for all these observables, including contributions from 
dimension-six operators. Finally, we compare them with the SM tree-level predictions, to 
identify the corrections from dimension-six operators. 

4.3.2 Fermion pair production at LEP2 

The observables are the total cross-sections and forward-backward asymmetries for fermion 
pair production, measured at different center of mass energies. The matrix element for 
e + e' ->//(/ ^ e) is given by 

M = 1 , A2 2 1 T* a 2 2 ^')7 M \9v - 9ll 5 ) u(p)u{k) lt , [g J v - ^ 7 5 v(k') 

(p + p') 2 -m 2 z + iY z m z A& w s 2 w K ' p V v A J 

^ 2 v{p>)Yu{p)u{k)lAk')i (75) 



(p + p' 
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where p, p' are the momenta of the incoming e + e~, and k, k' are the momenta of the outgoing 
fermions. The cross-sections and forward-backward asymmetries can be calculated from Ai, 
and Eq. ( 1591) can be applied to obtain the corrections from the operators. 

For f — e, there are additional contributions from the t-channel diagrams. The matrix 
element is 

M = ( I ^ ^ I t j^mY{9v-9ll 5 )u(p)mi,{9v-9 e Al 5 )v(k') 

(p + p') 2 - mi, + iT z m z AcfySw 

177/1 1 u(*)t* {9v ~ 9ll 5 ) uipW)^ (g e v - ^ 7 5 ) v(k') 



(p - k) 2 - m% + iT z m z Ac 2 w s 2 w 



A.TXOL Attq 

+ j^^v(p'h fl u(p)u(k)^v(k') - _ u(k)Yu(p)v(p')jAk )- (76) 
4.3.3 W mass 

The W mass is measured both at the Tevatron and LEP2. For the tree- level expression of 
mw, we first solve mw*ijn 2 v ) = m 2 ^, which gives 

m 2 w = m 2 W0 + Il ww {m 2 w ). (77) 

Combining Eqs. (l4"3"j) . (lo"Tj) and (j53"|) with Eq. ( 1771) . we find that the correction to the W mass 



is 



6m 2 w = n ww (m 2 w ) + 2 2 IW(0) - 2 3 2 n Z z(m|) + 4^Y^|n; 7 (0). (78) 

The W width is measured at the Tevatron. The tree level expression is 

3amw 



w 



4s 2 



(79) 



The correction can be calculated using Eq. (|59|) . 
4.3.4 DIS and atomic parity violation 

These are experiments performed at q 2 ~ 0. These low energy observables are usually 
expressed in terms of the effective couplings g(? and g^, which depend on s 2 ^. For the tree- 
level expressions, we will also include the factor p*(0), which is 1 in the SM, and takes the 
value of Eq. (1631) in the presence of new operators: 

P*(0) = 1 + Sp(0) = 1 - \n zz (0) + -4-IW(0). (80) 

m z m w 



The correction to an observable X is then given by 

g y tree 



5X = CiXf m6 = — 4^ 
ds 2 w 



Q x^tree 

Ss 2 w (0) + dXth 



p,(0)=l 



8p(0). 

p*(o)=i 



The observables include: 
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The weak charges for Cs and Tl, measured in atomic parity violation experiments. The 
weak charge is given by 

Q W (Z, N) = -2 \{2Z + N)C lu + (Z + 2N)C ld ] , (82) 

where Z and N are the proton number and the neutron number of the atom. The 
tree-level expressions for C\ u and C\ d are 

C lu = 2p*(0)g e A g$, C ld = 2p^)g\g d v . (83) 

The weak charge of the electron, Qw{ e ) , measured in polarized M0ller scattering: 

Q w (e) = -2C 2e = -ApM9 e A 9 e v (84) 

The effective couplings gi and gn for z/-nucleon scattering, measured at NuTeV. These 
are defined as 

9l = 9f >eS + 9t S , ffh = g%B + s1U- (85) 

where g^ eSl gR e s, gfes an d 9ReS are the effective couplings between the Z boson and 
the up and down quarks. The tree-level expressions are 

9l, eS = P*(0) 9 -^, 9i eS = P*(0) 9l i^, (86) 
*SU = *(°)^4 9U = PM 9 -^A- (87) 



• The effective couplings gy e and g v £ for v-e scattering, measured at CHARM II. The 
expressions are 

9v=PM9v, 97 = P*(V)9 e A- (89) 

4.3.5 W pair production 

So far we have been using the approach of Peskin and Takeuchi to study the effects of 
new operators. However, this approach only applies to processes involving light fermions as 
external particles, and cannot be used to study W pair production. Due to the relatively low 
statistics in the measurements, the constraints from W pair production are weak compared 
to other electroweak observables. Therefore we will ignore all loop effects, and only focus on 
the effects of operators Owb and . 

Using Eqs. (122))-(I25| and Eq. §3J, we have 



s 2 



if f Ana \ 
2 \ y V2G F ml J 

Note that Eq. (|53|) only has to do with the indirect corrections, so it is still valid. 



(90) 
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The operator Owb changes the mixing of the W 3 and B bosons. We then define sw 
and Cw* as the new mixing angle, 

> 2 _ 
'A 2 



n v 2 

Sw* — S W0 — ( -^WB~T2 S W c Wi C W* 



1 - S 2 

1 b W*i 



(91) 



so that the W + W Z and W + W 7 vertices from the kinetic term — \W / l lll W Iflu have the same 

form as in the SM. Note that this definition of s\y* is different from the one in Eq. (H0|) . 

(3) 

which was only valid for light fermions. In this way, the operators Owb and have the 
following effects: 



The SM Zff vertices are modified to: 



C 



zff 



Sw* c W* 



V 2 Sw 



'A 2 c w 



1 + c WBT - 2 — zj-f t 6 p l - s z w , ii + a 



J WB 



(92) 



where P L = (1 ~7 5 )/2. 



The SM W + W Z and W + W 7 vertices are modified. The contribution comes from 
Owb'- 

2 2 
Owb -> -igCwB^cwA^W+W- + igC WB ^s w Z^W+W; . (93) 



The W mass is changed to: 
1 



in 



w 



rrir 



2 



Ana 



V2G F m 2 z 



A 2 



or V 2 S W C W 1 (3) V 2 

' 2CwB T 2 ^s^-2 C ^ A^d^sJ 



c 



w 



-w °w 



(94) 



Using these results we can write down the matrix element. The process has a t-channel 
contribution Ait and s-channel contributions Ai^ and Aiz, which come from photon and Z 
boson exchange. They are given by 

-,2 



Mt 



—i 



e" 1 



2 s 2 



(95) 



M y = -ie 2 v(p')j p u(p)-^ x 



2 

(g^(k' - k) p - g^(q + k'Y + g pp {q + kf) + C W b^— (g P T ~ 9 vp <f) 

A z sw 



■>w* 



i + a 



WB 



v 2 s w 
A 2 c w 



-Pl - 1 + O 



'WB 



f C W 

A 2 s w 



V 2 Sw 



>w* 



u(p) 



q 2 — m 2 z 



X 



{gTW - k)" - g up (g + k'Y + 9 PP (q + kf) - Cwb^— (g pp q u - g up q p ' 

A z c w 



where p, p' are the momenta of the incoming e + e~ , q = p + p', k, k! are the momenta of the 
outgoing W + W~, and e* Al , e* A2 are the polarization vectors of the W + W~. q = p + p'. The 
cross section can be calculated from the matrix element. The correction due to Owb and 
Of is obtained by taking the linear part in C WB and Cf . Higher order terms in C/A 2 are 
neglected. 
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4.4 Total x 2 

In the calculation of x 2 , we choose the MS scheme, with the renormalization scale /i 2 = m|, 
even for processes in which the characteristic scale is not m 2 z . This is done for ease of 
calculation. The error introduced by this choice is of higher order and negligible. 

We find that the contribution from the operator is suppressed by the bottom quark 
mass, as can be seen from Eq. (I134p . Therefore we neglect this operator. The contributions 
from operators O b w and 0\>b also have a factor of mj. However, their effects can still be 
large, because the expressions contain the function bo{ml,ml,q 2 ), and its derivative with 
respect to q 2 is inversely proportional to m 2 : 



^ i. i 2 2 2^ 
-^b (m b ,m b ,q ) 



6m 2 



(98) 



Therefore, we will consider 10 operators: 

Owb, , Ofa, 0^, O^t, 0</>&, O tw , O bW , O tB , O bB - (99) 

Using Eq. f fl8|) . x 2 can be written as a quadratic function of Cf. 

X 2 = xLn + (Ci - C i )M ij {C j - Cj). (100) 

Here Xmm ^ s the minimum \ 2 m the presence of the new operators. Cj corresponds to the 
best fit value for Cj. 

In our calculation, we used the following input parameters: 

a(m%) = 1/128.91, G F = 1.166364 x lO" 5 GeV" 2 , m z = 91.1876 GeV, 

m t = 172.9 GeV, m b = 4.79 GeV. (101) 

We find xLn = §0.04, the 

Xmin P er degree of freedom is 0.80, compared with the SM value 
0.84. The matrix Mij and the best fit values Cj are given in Appendix [Bj 

4.5 Global fit 

The one-sigma bounds on the operators are given by x 2 ~ X m in = 1- By diagonalizing the 
matrix Mjj, we find ten linear combinations of C that are statistically independent. Their 
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best fit values and one-sigma bounds are given by: 



/ -0.961 
-0.064 
+0.268 
-0.008 
-0.014 
+0.009 
-0.016 
-0.004 
-0.001 
\ -0.001 



-0.273 
+0.159 
-0.940 
+0.019 
-0.065 
-0.107 
+0.030 
+0.008 
+0.001 
+0.001 



+0.029 
-0.701 
-0.063 
-0.095 
-0.336 
-0.326 
+0.137 
+0.034 
-0.003 
+0.505 



-0.004 
-0.680 
-0.182 
-0.086 
+0.344 
+0.322 
-0.137 
-0.034 
+0.003 
-0.505 



+0.024 
-0.015 
+0.088 
-0.003 
-0.389 
-0.590 
+0.128 
+0.039 
+0.007 
-0.689 



x 



A 2 



+0.000 
+0.130 
+0.002 
-0.991 
+0.017 
-0.009 
-0.000 
+0.001 
-0.000 
-0.000 

CwB \ 



\ 



r (3) 
W? 

C<f,b 
Ctw 
Cbw 

CtB 
CbB 



J 



+0.012 
+0.001 
-0.022 
-0.019 
-0.768 
+0.623 
-0.003 
-0.094 
+0.025 
-0.108 

/ - 



V 



-0.000 
-0.000 
+0.002 
+0.001 
+0.057 
-0.065 
+0.138 
-0.745 
+0.646 
+0.014 

-0.0004 
-0.013 
+0.011 
+0.59 
-0.22 
-1.76 
-2.2 
-9.2 
+ 102.4 
-1.36e+3 



+0.015 
+0.001 
-0.005 
-0.000 
-0.135 
-0.194 
-0.935 
-0.244 
-0.090 
+0.054 

±0.0029 

±0.014 

±0.023 

±0.27 

±1.10 

±1.63 

±11.9 

±21.1 

±50.4 

±1.38e+ 
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001 \ 


+0 


000 


-0 


001 


-0 


000 


-0 


039 


+0 


014 


-0 


229 


+0 


610 


+0 


757 


+0 


009 J 


\ 





TeV-^102) 



3/ 



where the 10 x 10 matrix is orthogonal. We can see that in the first row and the third 
row, the first two components are much larger than the other components. This means that 
these two rows approximately correspond to constraints on the coefficients Cwb and C^\ or 
equivalently, the S and T parameters. Since we are interested in the other eight operators, 
we can assume that Cwb and C^ always take the values that minimize the x 2 . This doesn't 

mean that we fix Cwb and C^\ The values that minimize the total x 2 depend on the values 
of the other eight coefficients. In other words, we let these two coefficients freely float. In 
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this way, we find the following constraints on the eight operators: 



— fl 709 


— n 701 


— n nnn 


-0.094 


-0.087 


-0.002 


-0.342 


+0.349 


-0.398 


-0.326 


+0.323 


-0.591 


+0.137 


-0.137 


+0.128 


+0.034 


-0.034 


+0.039 


-0.003 


+0.003 


+0.007 


+0.505 


-0.505 


-0.689 



A 2 



+0.128 
-0.992 
+0.017 
-0.009 
-0.000 
+0.001 
-0.000 
-0.000 

/ <t \ 

C^b 
Ctw 
Cbw 

CtB 
CbB 



\ 



-0.003 
-0.019 
-0.761 
+0.632 
-0.003 
-0.094 
+0.025 
-0.108 

/ 



V 



+0.000 
+0.001 
+0.056 
-0.065 
+0.138 
-0.745 
+0.646 
+0.014 

-0.011 
+0.59 
-0.23 
-1.75 
-2.2 
-9.2 
+102.4 
-1.36e+3 



-0.000 
-0.001 
-0.136 
-0.191 
-0.935 
-0.244 
-0.090 
+0.054 



±0.014 
±0.27 
±1.10 
±1.62 
±11.9 
±21.1 
±50.4 
±1.38e+3 / 



-0.000 \ 
-0.000 
-0.039 
+0.015 
-0.229 
+0.610 
+0.757 
+0.009 ) 

\ 



TeV 



(103) 



This is the main result of this paper. We can see that the first row is approximately a 
constraint on {Of^ + 0^)/\/2, which corresponds to the left-handed Zbb coupling, while the 
second row corresponds to a constraint on O^, which is the right-handed Zbb coupling. These 
are the tightest bounds, since the contribution arises at tree-level. The other constraints are 
mainly from loop-level effects. The third row is approximately a bound on Ctw- This can 
be compared with the bound obtained from the U parameter, Eq. f|T5|) . The results are in 
close agreement. The remaining five rows yield bounds on linear combinations of operators. 



5 Conclusions 

In this paper, we have studied the effects of non-standard top quark couplings on precision 
electroweak measurements. The top quark plays a role as a virtual particle in these measure- 
ments. Our study is based on an effective field theory approach, which allows us to calculate 
the self-energies of the electroweak gauge bosons at loop-level. 

We have examined the effects of nine dimension-six operators that generate non-standard 
couplings between electroweak gauge bosons and the third generation quarks. These opera- 
tors mainly contribute through loop corrections to the gauge boson self-energies, but some of 
them also have a tree-level contribution to the Zbb couplings. The contribution of one of the 
operators is suppressed by m h and is hence negligible. We have also included the operators 

(3) 

Owb and 0\ , in order to deal with the divergences that appear in our calculation. 

We have calculated the total \ 2 an d performed a global fit including these ten operators. 
We allow Cwb and C^ to vary, and thus obtain bounds on eight dimension-six operators. 
The result is shown in Eq. f 1 1 3 j) . The two tightest bounds are from tree-level contributions 
to the Zbb couplings, g\ and g b R , and the other bounds are from loop-level contributions. 
The best bound from loop-level contribution constrains ^ to be of order 1 TeV -2 . 



19 



If the bound on some linear combination of Cj is too weak, this bound cannot be trusted, 
because the linear analysis is not applicable if the coefficients are large. The contribution 
from dimension-six operators relative to the SM contribution is of order g 2 /(4ir) 2 x -nr- 
This should be less than unity, thus 

C 

-^«6.1x 10 3 TeV -2 (104) 

Even the weakest bound from our results does not exceed this limit, therefore all bounds 
can be trusted. 

Using Eq. (jlOOp . one can put constraints on a subset of these operators. For example, 
the one-sigma bound on the coefficient d, assuming only one coefficient deviates from its 
best fit value, is given by (7, ± M^ 1 ^ 2 , where M„ is the diagonal element of the matrix M 
and is not summed over i. 

We can also consider only one coefficient to be nonzero at a time. In this case, we found 
the constraints on each individual coefficient are: 



A 2 

U <f>q 

A 2 


A 2 
A 2 


= 0.016 ± 0.021 TeV" 2 
= 2.0 ± 2.7 TeV" 2 


(105) 
(106) 




A 2 " 


= 1.8 ± 1.9 TeV" 2 


(107) 




C<t>b 

A 2 


= -0.16 ±0.10 TeV -2 


(108) 




Ctw 
A 2 


= -0.4 ± 1.2 TeV" 2 


(109) 




Cbw 
A 2 


= 11 ± 13 TeV" 2 


(110) 




CtB 

A 2 


= 4.8 ± 5.3 TeV -2 


(111) 




ChB 

A 2 


= 8 ± 19 TeV" 2 


(112) 



Some of these operators are already constrained from other measurements. The operator 
Otw modifies the top-quark branching ratio to zero-helicity W bosons [12]. Recently, the 
combination of CDF and DO measurements of the P^-boson helicity in top-quark decays 
reports a measurement of /o, the zero-helicity fraction 



fo = 0.685 ± 0.057[±0.035(stat.) ± 0.045(syst.)] (113) 
This yields the constraint on C t w'- 

= 0.03 ± 0.94 TeV" 2 (114) 

A 2 

We see that Eqs. (|T5|) . (I109p . and (I114p give similar constraints. 



20 



Constraints may also be gleaned from B physics. The operators — O^j and 
affect the branching ratio for B — > X^. The following constraints are obtained in Ref. [23J: 

-1.6±1.3TeV~ 2 (115) 

0.030 ± 0.026 TeV -2 . (116) 

There are also contributions from O t w and Obw, but these are ultraviolet divergent. Thus 
there must be a tree-level contribution from another dimension-six operator, which masks 
the contributions from Otw and Obw- Therefore one cannot obtain bounds on these two 
operators. 

In addition, constraints on the operator O t w can also be obtained from Bd s — Bds mixing 

c 

= -0.06 ± 0.79 TeV" 2 . (117) 

A 2 

This is comparable to the bounds given in Eqs. (fT5l) . fllQQj) . and (11141) . There is also a 
contribution from 0$ — that is divergent and requires a tree-level contribution, so this 
operator cannot be bounded from B^ s — Bd,s mixing. 
Finally, we summarize all these constraints in Table [2j 



Coefficients 


Electroweak data 


W helicity 


B^X sl 


B d , s ~ B dtS mixing 


(C% + C%)/A> 


0.016 ±0.021 








( r {z) r^A /A2 
[^^g -^h ) / 


2.0 ± 2.7 




-1.6 ±1.3 




C^t/A 2 


1.8 ±1.9 








C # /A 2 


-0.16 ±0.10 








C(f,cf,/ A 2 






0.030 ±0.026 




C tw /A 2 


-0.4 ± 1.2 


0.03 ±0.94 




-0.06 ±0.79 


CbW 1 'A 2 


11 ± 13 








C tB /A 2 


4.8 ±5.3 








C bB /A 2 


8 ±19 









Table 2: Bounds on operators, in units of TeV 



We can see that for most operators our analysis gives the best bounds available, because 
electroweak data is the only way to access these operators so far. For the other operators, 
our bounds are comparable to the best bounds obtained from colliders and B physics. 
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A Gauge boson self-energies 



Here we give Hxy 
• {3) 

liww 



Uxy for all 9 operators. 

\q 2 -\{rn 2 t +m 2 b )^E 



= -N t 



n 



■7Z 



,(1) 



• Op 



g 2 v 2 

4vr 2 A 2 



-q 2 b 2 (m 2 ,m 2 b ,q 2 ) + ^ (m 2 &i(m 2 , m 2 , g 2 ) + m 2 &i(m 2 , m 2 , g 2 )) (118) 



n zz = -tv c 



_ 1 ^ 

COS 2 6>VF 47T 2 A 2 

2 . 



l -{l-sm 2 6 w )q 2 -- A {m 2 + m 2 )^E 
2 ™ 2 ^-Q-Isin 2 ^ 



n 77 = 



-g 2 (Q-^in 2 ^jfe 2 K,m 2 ,g 
+ ^ { m t b o( m t, m t, <? 2 ) + m l b o{ml, ml, q 2 )) 



2 6 w J b 2 (m 2 b ,m 2 b ,q 2 



c9 cos^8tt 2 A 2 



12 1 

-£ - -& 2 (™ 2 , m?, g 2 ) - -b 2 (m 2 b , m 2 , g 2 ) 



(119) 
(120) 
q 2 (121) 



n^vy — 
Ilzz = N, 



g 2 1 v 2 r /1 2 1 2 1 2 . 2 \ 

cos 2 % 4tt 2 A 2 [ V 4 4 18 / 



(122) 



-9 

+ 



n 77 = 
n 7 z 



i - I sin 2 0^ & 2 (^ 2 , m 2 , g 2 ) - Q - ^ sin 2 & 2 (™ 2 , ™ 2 , g 2 ) 
i (m 2 b {m 2 ,m 2 ,q 2 ) - m 2 b b {m 2 b , 
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r mr /r 



(123) 
(124) 



- (125) 



11^^ — 

n zz 
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, r <? 2 1 ^ 2 

AT — 2 

c cos 2 6 W 4tt 2 A 2 



- ^ 2 sin2 ^ £ 



n 77 = 
n 7 z 



(1 O7 / 2 2 2\ 2-2 /17 / 2 2 2\ 

-m t b {m t ,m t ,q ) - -q sm 6 w b 2 (m t ,m t ,q ) 

9 W 1 V 2 (I ,/ 2 2 2\\ 2 

^12^A^U ^'"^T 
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(126) 

(127) 
(128) 
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n Z z 



= o 



= N r 



1 v 2 



cos 2 9 W Air 2 A 2 



- ( \m\ - ^q 2 sin 2 9 W ) E 



[ ^m 2 b b (m 2 b ,m 2 b ,q 2 ) - ^q 2 sin 2 9 w b 2 (m 2 b , m 2 , q 2 ) 



n 77 = 



n 7 z = 



2 sin^/ 1 t; 2 /l 2 2 2 , 2 



(130) 

(131) 
(132) 
(133) 



R\vw = -N c g 2 j^—m t m b (E-b (m 2 ,m 2 b ,q 2 )) (134) 

Tl zz = (135) 

n 77 = (136) 

n 7Z = (137) 

IW = -N^^^-E-himlmlq^q 2 (138) 

Il zz = ~ N c9^^ (I - \ sin 2 9 w ^j (E - b (m 2 , m 2 , q 2 )) q 2 (139) 

n 77 = -N c g^ V ^sin 2 9 w (E-b (m 2 , m 2 ,q 2 ))q 2 (140) 
n 7 , = -"*^§£ (iH*^) (^<rnWM (14D 

IW = -N^^^E-bMlrnlq 2 )^ 2 (142) 

IW = -^^^Q-^sin 2 ^(E-6 K,m 2 ,g 2 ))g 2 (143) 

n 77 = -N c g^ V ^sin 2 9 w (E-b ( m lmlq 2 ))q 2 (144) 
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• o tB 



n 



■yZ 



• o bB 



Tlww 
U zz 



n 



77 



n 



7Z 



(146) 

AT y/2 vm t sm9w ( 1 4 . , \ , 2 2 2 , 2 

_JV ^A^3 Sin ^ C ° S ^ ^ ~ 6o(m * ' m * ' ^ ^ (148) 

-iVc^^ Q - I sin 2 ^) (E - m 2 , g 2 )) g 2 (149) 



(150) 



vm b sin 6V / 1 2 



4. 2 A 2 cos^U-i Sm2 ^) (^oK,^W (151) 



^c^^ sin ^ CO s^ ( E - ^oK,m 2 ,g 2 )) g 2 



N c g 



4tt 2 A 2 3 

\/2 vm b (\ 2. 2 

sin 6V 



4tt 2 A 2 V4 3 



(E-boi 



2 2 2\\ 2 

m b ,m b ,q )) q 



(152) 
(153) 



Here 0^ is the weak angle, N c = 3 is the number of colors. E = ^ — 7 + ln47r, and the 
functions bi are given by 



7 / 2 2 2\ 

b {m 1 ,m 2 ,q ) 
bi(ml,ml,q 2 ) 
b 2 (m\,m 2 2 ,q 2 ) 



1 , (1 — x)m? + rem 2 , — x(l — x)(/ 2 , 
In — f — dx, 

V 

1 , (1 — x)m 2 + xrrin — x(l — x)q 2 , 
x In — ^ —dx, 

1 , . , (1 — x)m 2 + xm 2 , — x(l — x)q 2 , 
x(l - x) In ^ '—kZ. i ^-dx, 

^ 



where /i is the 't Hooft mass. They have the following analytical expressions: 



bQ{m x ,m 2 ,q) = -2 + log — + 



mim 2 mi — m 2 



log 



mi 
m 2 



(154) 
(155) 
(156) 



+\y/\(rn 1 + m 2 ) 2 -g 2 ||(mi - m 2 ) 2 - q 2 \f(m\, m 2 , g 2 ), (157) 



where 



ill 2 2\ 



log 



-y/ (mi+tn 2 ) 2 -9 2 -i/ (mi-m2) 2 -? 2 



(mi+m 2 ) 2 -9 2 + i y (mi-m 2 ) 2 -9 2 

2arctan,/gtoJ 



g 2 < (mi — m 2 ) 2 

(mi - m 2 ) 2 < g 2 < (mi + m 2 ) 2 , 



^/g 2 -(m 1 -m 2 ) 2 -^/(/ 2 -(mi +m 2 ) 2 



(158) 
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and 

61 (m?, ml, q 2 ) 
b 2 {m\,m\,q 2 ) 



1 m i 1 
log — - 1 

/I 2 J q z 
m\{2m\ — 2m\ — q 2 ) 



1 1 

18 + 6 



m; 



2^2 

1 
6 



log 



- 1 



1 m\ — m\ + q" 

2 q 2 



&o(mi,m£j.i$) 



m 2 m^mj — 2m 2 — q 2 ) m 2 



H 2 



(q 



2\2 



m\ + m\ + g 2 



log^ 



b (m 1 ,m 2 , 4)160) 



B The matrix Mij and the best fit values & 



The matrix My and the best fit values Q in Eq. (llOOp are given by 
(1 TeV) 4 



M ■■ 

CwB 



A 4 



x 1(T Z 

^(3) 



c, 



0i 



c, 



06 



a 



fir 



C, 



611' 



IB 



a 



6B 



of 

og» 
(1) 

cv 

0<f,b 
Otw 
Obw 

tB 

o bB 



+1.10e7 
+3.06e6 
-3.16e5 
+5.47e4 
-2.70e5 
-6.16e3 
-1.35e5 
+3.11e3 
-1.71e5 
-1.40e4 



and 



+3.06e6 
+ 1.07e6 
-1.40e5 
-1.03e4 
-9.49e4 
+9.46e3 
-3.39e4 
+4.04e2 
-4.77e4 
-3.85e3 



-3.16e5 
-1.40e5 
+2.58e5 
+2.40e5 
+1.28e4 
-4.55e4 
+3.99e3 
-4.49el 
+4.96e3 
+4.35e2 



+5.47e4 
-1.03e4 
+2.40e5 
+2.39e5 
+1.16e3 
-4.42e4 
-1.28e2 
+3.21e0 
-8.20e2 
-3.34el 



-2.70e5 
-9.49e4 
+1.28e4 
+1.16e3 
+8.49e3 
-9.17e2 
+2.98e3 
-3.34el 
+4.21e3 
+3.40e2 



-6.16e3 
+9.46e3 
-4.55e4 
-4.42e4 
-9.17e2 
+9.83e3 
+ 1.13e2 
-1.46el 
+9.24el 
+3.20e0 



-1.35e5 
-3.39e4 
+3.99e3 
-1.28e2 
+2.98e3 
+1.13e2 
+1.78e3 
-5.16el 
+2.11e3 
+1.76e2 



+3.11e3 
+4.04e2 
-4.49el 
+3.21e0 
-3.34el 
-1.46el 
-5.16el 
+2.49e0 
-4.89el 
-4.42e0 



-1.71e5 
-4.77e4 
+4.96e3 
-8.20e2 
+4.21e3 
+9.24el 
+2.11e3 
-4.89el 
+2.67e3 
+2.19e2 



-1.40e4 
-3.85e3 
+4.35e2 
-3.34el 
+3.40e2 
+3.20e0 
+1.76e2 
-4.42e0 
+2.19e2 
+ 1.82el / 
(161) 



Ci 


C\YB 








C^t 








CtB 


CbB 


a/A 2 


+0.93 


-1.63 


-689 


+689 


+939 


-0.60 


+149 


+53.9 


-78.3 


+60.3 



(162) 



in units of TeV -2 . 

The numerical values of Ci depend on both the experimental values and the SM predic- 
tions. The matrix M is symmetric and positive definite, and its value only depends on the 
errors of different measurements. If any of the SM input parameters changes, the best values 
C{ will be affected, but the matrix M will not. The sizes of the one-sigma bounds on the 
operators only depend on the matrix M. 
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